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Abstract. The development of prebiotic homochirality on early-Earth or another 
planetary platform may be viewed as a critical phenomenon. It is shown, in the 
context of spatio-temporal polymerization reaction networks, that environmental 
effects — be them temperature surges or other external disruptions — may destroy 
any net chirality previously produced. In order to understand the emergence of 
prebiotic homochirality it is important to model the coupling of polymerization 
reaction networks to different planetary environments. 
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1. Introduction 

According to the prebiotic soup hypothesis (Wills and Bada, 2000), the 
early Earth had the supply of organic compounds needed to jump-start 
polymerization reactions that, through gradual complexification, led to 
the first biochemical networks displaying some of the characteristics 
attributed to life, such as metabolic activity and replication. Although 
the road is still obscure (Orgel, 1998), the situation is not all bleak. 
In 1953, Stanley Miller simulated tentative conditions of early-Earth in 
the laboratory to obtain amino acids from simple chemical compounds. 
That same year, Frank (Frank, 1953) proposed that auto-catalytic 
polymerization could explain the emergence of biomolecular homochi- 
rality, a clear signature of terrestrial and, possibly, all life (Bonner, 
1996): terrestrial amino acids belonging to proteins are overwhelmingly 
left-handed, while sugars are right-handed. 

If a bottom-up approach to the early development of life is adopted, 
the homochirality of life's biochemistry must have emerged dynami- 
cally, as reactions among the simplest molecular building blocks oc- 
curred with high enough yield. Alternatively, one may assume that, 
somehow, only monomers of a single chirality were present in the prebi- 
otic soup: they were made that way or brought here during the intense 
bombardment of Earth's infancy, that lasted to about 3.8 Gyr ago 
(Gomes et al, 2005). We would, however, still need to understand how 
homochirality developed elsewhere in the cosmos and not here, and 
whether it developed with the same chiral bias in more than one place. 
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Here, we consider the homochirality of life as an emergent process 
that took place on early-Earth and, possibly, other planetary plat- 
forms. As a starting point, we use the reaction-network model pro- 
posed by Sandars (Sandars, 2003), which includes enantiometric cross- 
inhibition. As shown in the interesting work of Brandenburg and Mul- 
tamaki (Brandenburg and Multamaki, 2004) (BM), Sandar's polymer- 
ization reaction network can be reduced to an effective spatio-temporal 
mean-field model, where the order parameter is the chiral asymmetry 
between left and right-handed polymers. To this, we add the effects of 
an external environment, showing that they can be crucial in the final 
determination of the net value of enantiometric excess, if any. 



2. Modeling Polymerization 

Sandars proposed the following polymerization reactions (Sandars, 2003): 
L n + L\ — ^ L n+ i; L n + R\ — L n R\; L\ + L n R\ — ^> L n+ \R\; 

ki 

and, R\ + L n R\ — ► R\L n R\, where k$ (ki) are the reaction rates 
for adding monomers of the same (opposite) chirality to a given chain. 
The network is completed by adding the four opposite reactions, that is, 
by interchanging L <-> R, and by adding a substrate S from which both 

left and right-handed monomers emerge, S k ^—> L\ and S k ^—> R l} 
where Cum determine the enzymatic enhancement of left and right- 
handed monomers, which is not known. We follow Sandars (2003) and 
choose Cl = [Ln], Cr = [Rn]- As remarked in BM, it is possible to 
truncate the system to N = 2 and still maintain the essential aspects 
of the dynamics leading to homochiralization. This allows us to model 
the reaction network as a mean-field theory exhibiting spontaneous 
chiral symmetry breaking (Weinberg, 1996). Our approach blends the 
work of BM with the pioneering work of Kondepudi and Nelson, where 
the reaction network was coupled to time-dependent external effects 
(Kondepudi and Nelson, 1985): chirality evolves spatio-temporally in 
contact with an environment. 

The equations can be simplified by assuming that the rate of change 
of [L2] and [R2] is much slower than that of [Li] and [R\] - The same for 
the substrate [S], so that d[S]/dt = Q — {Ql + Qr) — 0, where Ql and 
Qr are the source terms for monomers generated from the substrate 
[S\: Ql = k c [S](pC L + qC R ), and Q R = k c [S]( P C R + qC L ) (Haken, 
1983). The constants p = |(1 + /) and q = ^(1 — /) are given in 
terms of the fidelity of enzymatic reactions /, an adjustable parameter. 
As demonstrated by Kondepudi and Nelson (Kondepudi and Nelson, 
1983) and many others (Sandars, 2003; Wattis and Coveney, 2005; Saito 



capritc.tex; 5/02/2008; 9:44; p. 2 



Prcbiotic Homochirality as a Critical Phenomenon 



3 



and Hyuga, 2004; Gayathri and Rao, 2005; Avetisov and Goldanskii, 
1993), auto-catalytic networks show a tendency to bifurcate toward 
homochirality for values above critical, f c . The specific value of f c is 
model-dependent. 

Under the above assumptions, and introducing the dimensionless 
symmetric and asymmetric variables, S = X + Y and A = X — Y, 
where X = [Li](2fc 5 /Q) 1/2 and Y = [Ri](2k s /Q) 1/2 BM have shown 
that the polymerization equations reduce to 

dS_ 
~dt 

where Ao = (2ksQ) 1 ^ 2 , with dimension of inverse time. S = 1 is a fixed 
point: the system will tend towards this value at time-scales of order 
Ao- With 5 = 1, the equation for the chiral asymmetry has fixed points 
at A = 0, ± y/2f — 1, as pointed out in BM. An enantiometric excess 
is only possible if / > f c = 1/2. 



A c 7 1 ^ = l-^ (1) 



3. Introducing Environmental Effects 

We model the external environment via a stochastic spatio-temporal 
Langevin equation, rewriting the equations above as 

Ao_1 (f " W2S ) = i-^ + Ao^x,*) (3) 

A ° 1 (t£ " kv2ji ) = 2f s^ ~ SA + A °"^ (x ' t] ' (4) 

where k is the diffusion constant and £ represent white noise with zero 
mean and a two-point correlation function given by (£(x', t')£( x > t)) = 
a 2 (5(x' — x)5(t' — t), and a 2 measures the strength of the external 
influence. For example, in mean-field models of phase transitions, it 
is common to write a 2 = 2^ksT^ where fee is Boltzmann's constant, 
T is the temperature, and 7 is a viscosity coefficient. The equations 
can be made dimensionless by introducing to = Aot and xq = ^/Xo/kx, 
which determine the typical spatio-temporal scales in the system. The 
noise amplitude scales as Oq — > A ~ 1 (Ao/A:)' i / 2 a 2 , where d is the number 
of spatial dimensions. Using as nominal values for the parameters, 
ks ~ 10 _25 cm 3 s _1 , Q ~ 10 15 cm _3 s _1 , the diffusivity of water k = 
10 _9 m 2 s~ 1 , we obtain, t ~ (7 x 10 4 s)to and x ~ (lcm)xo. 

As in BM [see also (Saito and Hyuga, 2004; Gayathri and Rao, 
2005)], the concentrations are spatially-dependent quantities. This im- 
plicitly assumes that it is possible to define an effective correlation 
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volume within which the value of the chiral asymmetry A is fairly 
homogeneous. (We take S = 1.) Using well-known results from the 
mean- field theory of phase transitions (Landau and Lifshitz, 1980), we 
can easily compute the correlation length. From the equation for A, we 
obtain an effective potential V(A), 

V{A) = ^-f\n[A 2 + l] . (5) 

For / < 1/2, V(A) has a typical double- well shape, with minima at 
the fixed points A± = y/2f — 1. The correlation length, £, is given by 
£~ 2 (»4min) = V"(A m i n ), where A m m denotes a minimum of the poten- 
tial. For the fixed points, we get, £ 2 (U±) = f/(2f- 1). At f c = 1/2 the 
correlation length diverges, as it should for a critical point. However, 
the noise parameter a also controls the behavior of the system. Indeed, 
even if / = 1, an enantiometric excess may not develop if a is above a 
critical value a c . In analogy with ferromagnets, where above a critical 
temperature the net magnetization is zero, one may say that above a c 
the stochastic forcing due to the external environment overwhelms any 
local excess of L over R enantiometers within a domain of correlation 
volume V£ ~ £ d : racemization is achieved at large scales and chiral 
symmetry is restored throughout space. 



4. Numerical Results: Critical Point for Homochirality 

Salam (Salam, 1991) suggested that there should be a critical tem- 
perature T c above which any net homochirality is destroyed. However, 
he conceded that calculating T c would be quite challenging using the 
electroweak theory of particle physics. Here, we chose a different route 
which, we believe, will allow us to explore the qualitative aspects of the 
problem more effectively: the noise amplitude a may represent a sudden 
increase in temperature and/or pressure due to a meteoritic impact or 
volcanic eruption, or, possibly, due to a source of circularly-polarized 
ultraviolet light (Lucas et ai, 2005). 

The equation dictating the evolution of the enantiometric excess 
A, eq. 3, was solved with a finite-difference method in a 1024 2 grid 
and a 100 3 grid with St = 10~ 3 and 5x = 0.2, and periodic boundary 
conditions. In 2d, this corresponds to simulating a shallow pool with 
linear dimensions of I ~ 200cm. We prepared the system initially in 
a chirally pure phase, which we chose to be (A)(t = 0) = 1. The 
equation is then solved for different values of the external noise, a. 
As can be seen in Figure 1, for a 2 > a 2 ~ 0.65(£;/Ao) 3 / 2 , (A) — ► 0, 
that is, the system becomes racemized. (A) approaches a constant for 
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Figure 1. Average enantiometric excess versus "temperature" in three dimensions. 



large times, indicating that the reaction network reaches equilibrium 
with the environment. The results are ensemble averaged. For d = 2, 
a 2 ~ 1.15(fc/Ao). 

We can describe the environmental impact on homochirality at the 
microscopic level by introducing the "Ginzburg criterion" familiar of 
phase transitions (Landau and Lifshitz, 1980). Consider a correlation 
volume with (A) = 1 (or —1). What is the energy barrier (Eg) to 
flip half the molecules in the volume so that (A) — * 0? If is the 
number of molecules in a correlation volume, Eq = (N^/2)Ef , where 
Ef is the energy to flip one molecule. The Ginzburg criterion says that 
this energy is Eg — V^AV, where AV = \V(0) — V(±l)|. Comparing 
the two expressions we obtain, Ef = 2AV(V^/N^). From Equation 
5, |F(0) - V(±l)\ = 0.193. Now, ~ 4£ 3 = 4(fc/A ) 3/2 . (We set 
/ = 1.) Using for the microscopic spatial scale £ m i cro — (Q/ks)^ 1 ^ 6 , 
and that N% ~ (^/^micro) 3 , we obtain Ef ~ 1.5 x 10~ 26 m 3 . [The energy 
has dimensions of (length) 6 *.] To complete the argument, we use that 
the critical "environmental" energy to restore the chiral symmetry was 
obtained numerically to be (cf. Figure 1), E env ~ 0.65(/c/Ao) 3//2 = 
0.65 x 10" 6 m 3 . We thus obtain the ratio, £//£ e nv| ^ 2.3 x 10" 20 . 

It is tempting to compare this result with possible sources of ho- 
mochirality. For example, weak neutral currents are expected to induce 
an excess at room temperature of g = AE/ksT ~ 10~ 17 (Kondepudi 
and Nelson, 1983). In the language of the present work, they would 
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induce a tilt in the potential V(A) proportional to g. Thus, within the 
violent environment of prebiotic Earth, effects from such sources, even 
if cumulative, would be negligible: any accumulated excess could be 
easily wiped out by an external disturbance. Further work along these 
lines is in progress. 
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